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Influence of disorder strength on phase field models of interfacial growth 
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We study the influence of disorder strength on the interface roughening process in a phase- 
field model with locally conserved dynamics. We consider two cases where the mobility coefficient 
multiplying the locally conserved current is either constant throughout the system (the two-sided 
model) or becomes zero in the phase into which the interface advances (one-sided model). In the 
limit of weak disorder, both models are completely equivalent and can reproduce the physical process 
of a fluid diffusively invading a porous media, where super-rough scaling of the interface fluctuations 
occurs. On the other hand, increasing disorder causes the scaling properties to change to intrinsic 
anomalous scaling. In the limit of strong disorder this behavior prevails for the one-sided model, 
whereas for the two-sided case, nucleation of domains in front of the invading front are observed. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Interface growth in disordered systems has been a sub- 
ject of great interest in non-equilibrium statistical physics 
for more than a decade. Many different examples of inter- 
faces growing in heterogeneous systems have been found 
in phase separation by nucleation and growth [l|,|2j, solid- 
ification [J Q , and fluid invasion in porous media [E Q , 
among others. 

Phase field models of increasing complexity have in 
recent years been extensively used in studying inter- 
face roughenin g a s well as microstructure formation 
[E IE B M EETI1 El El- A particularly interesting 
problem of interface roughening is that associated with 
a fluid invasion front moving into a disorder medium Q , 
which can be experimentally studied with the Hele-Shaw 
cell set-up [II, EE EE EE El- In modeling such a fluid in- 
vasion experiment, two different ways to consider the mo- 
bility parameter in a Model B type of phase field model, 
called the one-sided and symmetric models, were use d by 
Hernandez-Machado et al. [lj| and Dube et al. [20, |2l[ , 
respectively. The difference between these two cases is 
that in the two-sided model the mobility is constant 
throughout the system, while in the one-sided model it 
becomes zero in the phase which is being invaded. 

In this paper our aim is to carry out a detailed anal- 
ysis about the influence of the strength of the disorder 
in the two types of models described above. Both cases 
can be described by a generalized Cahn-Hilliard equation 
or Model B with quenched disorder in the background 
medium. The boundary conditions are used to couple 
the system to a reservoir of the invading phase with a 
constant mass flux. At the linear level of small front 
fluctuations, both phase field models can be analyzed 
through linearized interface equations, that is, the evo- 
lution of the interface can be described in terms of the 
interface profile alone. Moreover, the bulk diffusion fields 



implicit in this description cause the interface equation to 
become spatially non-local [IE HE HH . It is thus of inter- 
est to examine how the models are influenced by varying 
disorder strength, which can be easily realized in the ex- 
periments, too. To this end, we first define a critical value 
of the disorder strength cr c , above which the disorder be- 
comes strong. We find that at weak disorder strengths 
(<r < <7 C ) both models have the same interface scaling 
behavior given by super-rough anomalous scaling, which 
is also predicted by the linear theory. We observe clear 
deviation from this scaling when the disorder strength 
comes close to cr c , where super-roughness disappears and 
weak intrinsic anomalous scaling arises. Furthermore, in 
the limit of strong disorder (a > a c ), the two models are 
no longer equivalent. The one-sided model can still be 
applied to describe diffusive liquid invasion in good agree- 
ment with experimental results of Ref. [16j |. However, in 
this limit the symmetric model exhibits nucleation of the 
invading phase in front of the advancing interface. 

The structure of the paper is as follows. In Section |TT] 
we introduce the two versions of the phase field model, 
consider the linearized interface equation (LIE) valid in 
the weak disorder limit, and give an estimate for the 
strong disorder limit, thus introducing a scale for the dis- 
order strength. Section IIIII defines the concepts of scal- 
ing in interface roughening, including super-rough and 
intrinsic anomalous scaling. In Setion HVl we present our 
numerical results, and finally give our conclusions in Sec- 
tion IVl 



II. THE PHASE FIELD MODEL 

The phase field model we are using describes a system 
of two phases separated by an interface. We introduce 
a locally conserved field <f> to represent the two phases 
of the problem with the equilibrium values <f> e = ±1 in 
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such a way that the interface position is at 4>{x, h) = 0. 
The dynamics of the field is then assumed to follow a 
conserved equation based on a time-dependent Ginzburg- 
Landau Hamiltonian (model B Q): 

dc/>/dt = -V-j, (1) 

where the current is proportional to the gradient of the 
chemical potential j = — A/(0)V/x. Here, the chemical 
potential is /i = SJ-'/Scj) and the free energy is taken to 
be of the form T[<j>] = f dr(V((f>) + (eV0) 2 /2). A sim- 
ple double-well potential is chosen with a linear random 
term: 

V{4>) = -~f + -r?(r)4> (2) 

The variable rj(r) is taken to be stochastic and it plays 
the role of spatially quenched disorder in the system. 
The disorder is characterized by its average (77), and its 
standard deviation a, which characterizes the disorder 
strength. The disorder also has a correlation length l corr , 
which in a numerical scheme is most conveniently set to 
the lattice spacing. Note that by considering 77 as local 
average over area of size Z corr , the choice of the length 
Icorr also enters the disorder strength. That is, the stan- 
dard deviation of 77 is a when observed at scale l CO rr- 
The surface tension of this model can be calculated with 
the disorder-free "kink" solution <pQ, or Goldstone mode, 
given by the lowest energy solution of the boundary con- 
ditions <f>{— 00) = —1 and </>(+oo) = 1. The resulting 
dimensionless surface tension is 7 = V2/3 ~ 0.47 @. 
The equation for the dynamics of the phase field reads 

^ = VM (0)V/z 

= VM(0) V [ - <p + 4> 3 - e 2 V 2 - t}(t)] , (3) 

where M(cj>) is a mobility parameter which may depend 
on the phase field. In the sharp interface limit e — > 0, 
the normal velocity of the interface can be obtained by 
integrating equation ([3]) in a region around the interface, 

v n ~ j n \ ± = M((f>)d n ^\ + - M(<j>)d nf i\-, (4) 

where the subscripts + and — correspond to the two 
phases of the system. In our study, two different func- 
tional forms for the mobility M{4>) will be considered. 
First, we assume a symmetric parameter M = Mq con- 
stant for the whole system independent of the field </>. 
In this case, the velocity of the interface is controlled 
by the difference between incoming and outgoing cur- 
rents. Second, we will consider a one-sided parame- 
ter M = Mq6 ((/>), 9(cf>) being the Heaviside step func- 
tion, which is zero in the phase where <j) < 0. Note 
that the normal velocity of the interface then reduces 
to v n ~ — Mo<9 n /i|_, that is, it is only proportional to the 
outgoing current. As we will see in our numerical results, 
the two models can give different results and describe dif- 
ferent physical situations depending on the strength of 
the disorder. 



Here we consider the so-called driven case, where the 
mean velocity of the interface is fixed to a constant value. 
The relevant boundary condition is to impose a fixed gra- 
dient of the chemical potential at the bottom of the sys- 
tem, V/i = — Vy. Combined with the initial condition 
of 

/ +1 y < H init ; , , 

<t> = \-ly>H inU , ^ 

the boundary condition leads to phase <fi = +1 invading 
phase <p — —\ and the interface moving with a constant 
average velocity. In the general context of the phase field 
model we will refer to these phases as A and B, respec- 
tively. In terms of liquid front invasion into a Hele-Shaw 
cell, these phases would be liquid and air, respectively, 
whereas in terms of phase separation they would be the 
phases rich in components A and B. 



A. The linearized interface equation 

The dynamics of a front in the phase field model de- 
scribed above can also be considered in terms of an 
integro-differential equation for the interface, which re- 
duces the dimensionality of the problem by one, obvi- 
ously a desirable property. However, this description can 
only be obtained as a perturbation expansion around a 
flat front with small fluctuations, and since the fluctua- 
tions are caused by the disorder, this also means weak 
disorder. Because the front dynamics are only obtained 
to first order in small fluctuations we call it the linearized 
interface equation (LIE). It is noteworthy that the long- 
ranged effects of mass conservation in the phase field re- 
sults in the LIE being spatially non-local, even though it 
is linear. This means that the LIE takes the form of con- 
volutions in position, which can be made local in Fourier 
space. 

The LIE can be obtained using the Green's function G 
to express the phase field equation, Eq. ([3]), in terms 
of an integro-differential equation [20j, and then pro- 
jecting it to an interfacial description with the operator 
= / dud u (f>o(u)[-}. The LIE takes the form of dis- 
persion relation for Fourier components of small front 
fluctuations h(k, t) around the average interface height 
H (t), {i.e. H(x, t) = H {t) + h(x, t)): 

d t h(k, t) = -(H \k\ + 7 |£f t) + \k\fi(k, t), (6) 

where 7 is surface tension, Hq — V is the (constant) av- 
erage front propagation velocity, and the disorder term 
is the Fourier transform of the two-dimensional (2D) dis- 
order along the front, or fj{k 1 1) = J dxe~ lkx rj(x, H(x, t)). 
From the dispersion relation ^ , one immediately obtains 
the crossover length scale 
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when the two dispersion terms are equally strong. Phys- 
ically the dominant dispersion mechanism then changes 
from surface tension to mass transport. In addition, be- 
cause the conservative character of the capillary disorder 
term \k\r)(k,t), it turns out that the crossover length £ x 
acts as an upper cutoff for the correlation length of fluc- 
tuations. This means that in the long wavelength region, 
where mass transport controls the dissipation of front 
fluctuations, the interface is asymptoticallyflat. This ef- 
fect has been numerically shown in Refs. [5|, [2(| l2ll I22I . 
and also in a general context of kinetic roughening [23j . 



B. Definition of disorder strength 

In order to study the influence of disorder strength in 
the phase-field model of Eq. (3), we need to define a 
measure for the relative strength of the disorder. This 
can be achieved by comparing the disorder contribution 
in the dimensionlcss bulk free energy of Eq ([2]) to the 
surface energy. We do this by considering a domain of 
linear size r, where the local disorder average (r}) r is a 
stochastic variable with standard deviation a r oc r a, 
where a is the standard deviation of a single disorder 
site, which is of linear size l CO rr- The underlying disorder 
then has a correlation length l corr . Considering a circular 
domain of radius r, it is energetically favorable for this 
domain to be of the opposite phase than its surroundings 
if 



(8) 



where the miscibility gap in our model is A(f> — 2. The 
LHS is the energy cost of the interface, whereas the RHS 
is the energy gain due to disorder. We consider the local 
disorder average (r)) r on the fluctuation site to be as large 
as its standard deviation a r — a l CO rr/ \f^r. Then Eq. ([5]) 
gives the condition 



v > -j 7, 



(9) 



for when the disorder can locally dominate the bulk en- 
ergy and thus is defined to be strong. The order of mag- 
nitude estimate for strong disorder in our dimensionless 
units {lcorr = 1, e = 1) is thus obtained as the variance 
being of the same order as the surface tension a — a c ~ 7. 
Note in particular that no r dependence remains in the 
estimate, and thus the relative disorder strength will be 
the same at all length scales (larger than the interface 
width e and disorder site size l cor r)- 



III. SCALING OF ROUGH INTERFACES 

The statistical treatment of a ID interface H(x,t) is 
usually done by studying the scaling properties of its 
fluctuations over the whole system of total size L [24j. 
For scale-invariant growth the lateral correlation length 



of the fluctuations is expected to grow in time follow- 
ing a power law £ c ~ t x l z until it reaches the sys- 
tem size L, defining a saturation time t s ~ L z . Al- 
ternatively, the global width of the interface W(L,t) — 

{[Hix^-H] 2 ) 1 ' 2 increases as W(L,t) ~ t p for t < t s 
and becomes constant W(L,t) ~ L x for t ^ t s . Here, (..) 
denotes average over different noise realizations and the 
overbar is an spatial average in the x direction. The 
quantities x, and z are the roughness, growth and 
dynamical exponent, respectively, and they are related 
through the scaling relation x = 0Z. In the standard 
Family- Vicsek scaling framework [25[ , this set of scaling 
exponents fully describe the scaling properties of the in- 
terface fluctuations. 

However, experimental results and several growth 
models have appeared in the last decade showing that 
global and local scales are not always equivalent. This 
is called anomalous scaling (26l. |27|. Therefore, one has 
to compute the interface width at different window sizes, 
w(£,t) = ({[H(x,t) - (iJ)f] 2 )^) 1 / 2 , where (,..) e denotes 
an average over x in windows of size £ < L. For scale- 
invariant interfaces local fluctuations are expected to in- 
crease as 



w(£,t)=£*g(£/t 1 /*), 

with the corresponding scaling function 
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u -(x-xioo) f or u K< 1; 
vT x for u >> 1, 



(10) 



(11) 



where xioc is the local rough exponent and it character- 
izes the roughness at small scales. One of the implications 
of anomalous scaling is that the local width saturates 
when the correlation length reaches the system size, i.e. 
at the time t s and not at the local time tg ~ £ z as occurs 
in the Family- Vicsek scaling. There is an intermediate 
regime between te and t s where the local width grows as 
w(£, t) ~ t p * with 0* = p - xioe/z. 

Another useful technique to determine the whole set 
of scaling exponents is to study the power spectrum of 
the interface S(k,t) = (h(k,t)h(—k,t)}. In the presence 
of anomalous scaling it is expected to show the following 
scaling 



S{k,t) = k- {2x+1) s A {kt 1/z ), 
where the scaling function has the general form 



sa(u) 
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Hx-xs) 



for 
for 



u » 1; 
u « 1, 



(12) 



(13) 



and Xs defines the spectral roughness exponent. Different 
scaling arises from this scaling function [26[. For Xs < 1 
it is always valid that xioc = Xs and the Family- Vicsek 
scaling is recovered whenever xioc = X- I n contrast, the 
so-called intrinsic anomalous scaling is observed when 
Xs 7^ X- Note that for this kind of anomalous scaling 
a temporal shift in the power spectrum is observed. The 
quantification of this shift based on the above scaling 
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FIG. 1: Interface widths for the two models at different disorder strengths (in dimensionless units). Results from the symmetric 
and one-sided model are given in the upper(a,b,c) and lower(d,e,f) panels, respectively. Disorder strength is varied from left to 
right as: weak disorder (a — 0.2), intermediate disorder (a = 0.5), and strong disorder [a = 1). Fitted growth exponents for 
the smallest and largest slopes are given in the figures, with solid lines representing the fits. 



form is unreliable and inaccurate, however, and thus we 
refrain from giving a measure for x~ Xs m the one case 
we observe scaling of this type. On the other hand, for 
Xs > 1 Xioc — 1 and the super-rough scaling appears 
when x = Xs- 

It is worth to remember here that the crossover length 
£ x , present in the imbibition phenomenon becomes a very 
important scale in the kinetic roughening process. More 
precisely, extensive numerical studies [1, Ho, [HI, [22j have 
shown that fluctuations of the interface do not evolve 
in time beyond this crossover length. Therefore, the in- 
terface fluctuations reach the steady-state at t s — £J 
instead of saturating at L z . Additionally, since the in- 
terface is asymptotically flat at scales larger than the 
crossover length, the crossover shows up as a maximum 
at the corresponding wavevector in the saturated struc- 
ture factor S(k,t — * oo). 

IV. NUMERICAL RESULTS 

Our study will be focused on the influence of the disor- 
der on the scaling behavior of the fluctuating interface for 
the two cases of the mobility M in the phase field model, 
that is, the one-sided model and the symmetric case. In 
our simulations we have used gaussian distributed dis- 
order with zero mean (r/) = 0, and different disorder 
strengths (standard deviation of r/) a. In the numer- 
ical scheme the disorder will have a correlation length 
as long as the lattice spacing, meaning that the lattice 



spacing normalizes the standard deviation when dimen- 
sionless numbers in the numerical scheme are turned to 
physical units. 

For the driven boundary condition, the methods used 
to obtain the interface description, and ultimately the 
LIE, predict that the mean value of the disorder will be 
irrelevant, as it has no contribution for the interface prop- 
agation. This has been numerically verified by our simu- 
lations of the full phase field models using different values 
for the disorder average. Here we only report results with 
(r?)=0. 

As numerical method, we chose the simple explicit Eu- 
ler scheme, where the disorder can be straightforwardly 
added. In two dimensions the timestep requirement 
of this method is not too restrictive for systems large 
enough for present consideration. 

A. Weak disorder 

For weak disorder, both models of the phase field (one- 
sided and symmetric) are expected to be equivalent, since 
the same LIE describes both cases in this limit. The weak 
disorder regime corresponds to d C 7, where 7 is the 
surface tension. In our dimensionless units (e = 1 in Eq. 
©) 7 ~ 0.47. Numerically the roughness and growth ex- 
ponents of Xs — 1-3 and (3 ~ 0.4, within the super-rough 
scaling description (xioc — L X — Xs), were observed 
by integrating the LIE 22]. This is in agreement with 
our results for both phase field models at weak disorder, 
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which are shown in Figs. Q] and El on the left. From 
Fig. Ufa) and (d) we observe x — Xsj since no temporal 
shift in the structure factor is present (unlike Fig. Etc)). 
Moreover, in Fig. [3J we observe that at small disorder 
strengths a < 0.2 in accordance with a -C 7, the spectral 
roughness exponent saturates to the LIE value, and is 
independent of the disorder strength. 

We also find that numerical artifacts in the interpre- 
tation of the interface from the phase field model appear 
when the disorder is so small that the global interface 
width W(L,t) is much less than the numerical lattice 
constant. This is due to the interpolation required to 
obtain the location of the interface between two sites 
of the numerical grid, and shows up as prominent pe- 
riodic oscillation in interface width with oscillation time 
Ax /Hq- While this numerical artefact can be removed 
by decreasing Ax, it means that in the convenient and 
typically used 0, [U, [22j dimensionless units, which 
lead to Eq. ([3]) with e = 1, only a very limited disorder 
strength range leads to the universal weak disorder limit. 
This range is roughly 0.1 < a < 0.2, when Ax = 1. 



B. Strong disorder 

In the regime of strong disorder, a > 0.5, different 
scenarios occur depending on which model is used. Be- 
low, we discuss the cases of the symmetric and two-sided 
models separately. 

1. Symmetric model 

According to the analysis of the disorder strength 
above, in this limit it is favorable for the phase field model 
to spontaneously create disperse domains of one phase 

(A) within the region that is initially of the other phase 

(B) . Droplets of phase A will form in phase B, and this 
mixture will initially cover most of the system. However, 
local mass conservation must still be valid regardless of 
any nucleation events. This means that mass must be 
diffusively transported from the A phase to the location 
where it nucleates within the B phase to facilitate a grow- 
ing droplet. 

This is exactly what happens in the symmetric model 
(see Fig. 2]), when one imposes the initial condition of Eq. 
([5]). Since there is no characteristic scale for the domain 
creation, the droplets are not restricted by the crossover 
length £ x , which acts as a cutoff for the interface fluctua- 
tions and therefore, the surface roughening at large scales 
is different respect to a weaker disorder strength. This is 
observed in Fig. Etc), where the interface power spectrum 
is plotted for a higher disorder strength at different times 
using the symmetric model. We can see that the fluctu- 
ations are not saturated, indicating that the crossover 
length £ x (represented by the dashed line) is not acting 
anymore as an upper cutoff. In addition, our numerical 
results for the symmetric model at strong disorder show 



three differences to the weak disorder case. First, the lo- 
cal growth exponent (3* approaches the global exponent 
(3 (see Fig. djc)). Second, the spectral roughness expo- 
nent decreases drastically to the range of \s — 0.5 (see 
Fig. [3j) and third, a temporal shift appears in the power 
spectrum (see Fig. Etc)). We can thus conclude that the 
scaling picture of interface fluctuations changes from su- 
perrough to intrinsic anomalous scaling, where Xs 7^ X- 

The picture becomes problematic for strong disorder, 
however, because the interface becomes less and less rep- 
resentable by a single valued function H (x, t). This is due 
to the interfacial area becoming more and more tattered 
by overhangs, droplets and bubbles. This also means that 
some numerical tricks are needed to distinguish the inter- 
face from these bubbles and droplets. This distinction is 
essentially made by finding a path for the phase bound- 
ary across the system that locally has as small height 
jumps as possible. This works relatively well as long as 
the disorder is not much stronger than a = 1 in our di- 
mensionless units. However, note anomalous fluctuations 
in Fig.[Ttc) around value 2.6 on the vertical axis, that are 
due to the abovementioned reason. 



2. One-sided model 

Using the one-sided model allows us to suppress the 
domain creation in phase B, where the mobility param- 
eter is zero. Then, the position of the interface H (x, t) 
can be found by taking the largest height where the phase 
A has advanced to at time t, coming from the phase B 
when the phase field is above zero. In Fig. Q] we show an 
example of the interface profile at different times for a 
strong disorder, a = 1.0. 

The growth exponent (3 ~ 0.5 measured for strong dis- 
order strength (see Fig. [Tte) and (f)) agrees with the 
experimental value of (3 — 0.50 ± 0.02 reported in Ref. 
[16( for liquid front dynamics into a Hele-Shaw cell. Like- 
wise, a similar variation in the spectral roughness expo- 
nent Xs, which changes from Xs — 1-23 to Xs — 0.91 when 
the disorder strength is increased (see Fig. [3j, was also 
experimentally observed in the same reference, whith a 
variation from Xs = 1-1 ± 0.1 to Xs = 0.9 ± 0.1 when 
the capillary forces of the Hele-Shaw cell were increased 
(see Fig. 15 in Ref. [16]). On the other hand, we numer- 
ically observe that the crossover length £ x still acts as a 
cutoff length for the interface fluctuations at strong dis- 
order (see Fig. E])- These results indicate that the model 
can still describe the imbibition phenomenon at strong 
disorder. 



V. CONCLUSIONS AND DISCUSSION 

In this work, we have studied two different ways of 
considering the influence of the mobility parameter in 
a Model B type of phase field model with a Ginzburg- 
Landau type free energy. The main experimental context 
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FIG. 2: Structure factors at ten equidistant time intervals for the two models at different disorder strengths (in dimensionless 
units). Results from the symmetric and one-sided model are given in the upper(a,b,c) and lower(d,e,f) panels, respectively. 
Disorder strengths are varied from left to right as weak (a = 0.2), intermediate (a = 0.5) and strong (a — 1). Fitted roughness 
exponents are given in the figures, with solid lines corresponding to the fits. The dashed vertical line corresponds to the 
crossover point k x = 27r/£ x as obtained from the LIE. 
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FIG. 3: Spectral roughness exponent \s for both phase field 
models as a function of disorder strength. 



considered here is liquid front invasion into a Hele-Shaw 
cell with quenched disorder [lq . We have focused on 
the case of driven front invasion, where there is a forced 
constant mass flux into the system that follows locally 
conserved dynamics. The symmetric model, studied for 
example in Refs. [f| HHHll], uses a constant mobility fac- 
tor, whereas the one-sided model, studied for example in 
Refs. [ill, [28j], uses a mobility that is zero in the receding 



phase, which we call phase B. 

We note that both models have previously been turned 
into non-local linear interface equations (LIEs) in the 
limit of small front fluctuations, which is equivalent to 
weak disorder [lj| [U [22j • These LIEs are identical for 
both models, and therefore both models are expected to 
have identical scaling behavior at the weak disorder limit. 
This is verified by direct comparison of numerical simu- 
lations. Furthermore, these results also agree with the 
relevant Hele-Shaw experiments fl6j |. 

We give an estimate for the strong disorder limit by 
comparing the disorder contribution to bulk energy to 
the surface tension. We find that the linear weak disor- 
der limit is found well below this disorder value, and that 
only in this limit does the roughness exponent not contin- 
uously depend on the disorder strength in either model. 
This means that a well-defined region of universality only 
exists at the weak disorder limit. 

Numerically we study the dependence of the growth 
and roughness exponents of the invasion fronts as a func- 
tion of the disorder strength. Our results are consistent 
with a (continuous) change of scaling behavior from su- 
perrough to intrinsic anomalous scaling, when the disor- 
der strength is increased from weak to strong. 

At strong disorder, the symmetric model is no longer 
found to correspond to the Hele-Shaw experiment due 
to domain creation of the invading phase in front of the 
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FIG. 4: An example of a set of rough fronts for strong dis- 
order cr = l using both models: one-sided model (top) and 
symmetric model (middle). The bottom figure shows points 
of zero <j> at fixed time in the symmetric model, demonstrating 
the nucleating domains. 



propagating interface. As our analysis shows, the do- 
main growth can occur at all length scales (larger than 
the disorder site size l CO rr of the system) without any 
characteristic radius, which is a phenomenon observed in 
other experimental situations such as nucleation on dis- 
locations [1^] or binary mixtures 0. In contrast, the re- 
sults for the one-sided model do agree well with the Hclc- 
Shaw experiments [161 ]. even as the disorder strength is 
increased. 

We hypothesize that the change in scaling behavior 
is due to decreased effect of surface tension and mass 
transport in interface roughening as the disorder becomes 
strong. Since the argument of strong disorder is the dom- 
inance of disorder over surface tension, the scaling in the 
regime dominated by the surface tension (characterized 
by having a correlation length l c < £ x ) should change for 
both models at strong disorder. In the one-sided model, 
mass transport still restricts interface roughening, and 
thus the crossover scale £ x persists. Conversely in the 
symmetric model the nucleated domains create roughen- 
ing by avalanches that are not controlled by mass trans- 
port from the reservoir, and thus the crossover scale £ x 
becomes irrelevant. 
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